Projected p-wave superconducting wave- functions for topological orders 
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In this paper we develop a sysmatical theory for topological orders by the projected p-wave 
superconducting (SC) wave-functions and unify the different topological orders for spin models into 
the fermionic picture. We found that the energy for the fermions at k = (0,0), (0,7r), (7r,0), (vr,7r) 
acts as a topological invariable to characterize 16 universal classes of different topological orders for 
the spin models with translation invariance. Based on the projected p-wave SC wave- functions the 
topological properties for the known topological orders in the exact solved spin models are obtained. 
Finally new types of topological orders are predicted. 
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Topological order is a new paradigm that lies be- 
yond the realm of Landau's theory based on the concept 
of symmetry-breaking for classical statistical systems. 
Physically, the topological orders describe the quantum 
orders with the long range entanglements in a gapped 
quantum ground state. A well-known example for the 
topological order is the fractional quantum Hall (FQH) 
state, which carries fractional charge and anyon statis- 
tics. To learn the nature for the FQH states, Laughlin 
wave function is a good staring point. Recently, many 
spin models with different topologically ordered states 
were found, such as the toric-code model [l|, the Wen's 
plaquette model [2, l3| and the Kitaev model on a hexag- 
onal lattice [J]. It is interesting that all the topologi- 
cal orders of the exact solved spin models can be de- 
scribed by p-wave superconducting (SC) wave functions 
[3|, [sj, la, [tI, y. People may ask two questions : if the 
topological orders for the different spin models arc the 
same one? and how to use p-wave SC wave functions 
to describe the topological properties of the spin mod- 
els? In this paper, we would like to develop a sysmatical 
theory for the topological ordered states based on the 
(projected) p-wave SC wave-functions and learn the na- 
ture of the different kinds of topological orders for the 
spin models. 

The general "mean-field" spinless fermion Hamilto- 
nian - We start with the general "mean-field" spinless 
fermion Hamiltonian with nonzero p-wave SC order pa- 
rameter 

Hjnea7i = ^ ^!"u V'j + '^{i'huji'] + h.C.) + A^ ^ V'lV'i- 
ij ij i 

(1) 

where Uij, rjij are 2 by 2 complex matrices, fi is the 

chemical potential. Let l^meaV?' } be the ground state 
of Hmean- To represent the topological order with 
full gapped excitations, each site must have just sin- 
gle fermion. Then we add a constraint to the p-wave 
SC wave functions. The topological ordered state can 
be obtained from the mean- field SC state l^meon' } by 



projection it into the subspace with single fermion per 
site, l^b'j''"') = P|*mea«''^). Here one has the projec- 



' spin 



tion operator as P = J|,- 



i-(-i) 



Because the fermion 



number can only he or \ on each site, the constraint 
for the fermion number on each site can be reduced into a 
constraint for total fermion numberl Because the fermion 
number for the superconducting state is not conserved, 
we can only distinguish the ground state with an even to- 
tal fermion number or an odd one. On the even-by-even 
(e*e), even- by-odd (e*o) and odd- by-even (o*e) lattices, 
the total fermion number N must be even, the projection 



operator turns into P = P^ = 



1+ 



On an odd-by- 



odd (o*o) lattice, the total fermion number must be odd, 
the projection operator P is reduced into P° = — 



Here N = 'J2''Pi''Pi = X) V'kV'k is the total fermion num- 

i k 

ber operator. 

For the translation invariant ansatz, one can rewrite 
the above "mean-field" fermion Hamiltonian in momen- 
tum space as H^ean = Z]k*k^(l*)*k with *k = 



^t ) a,nd *, 



( V'k V'-k ) • We note that 

ii^-ki V'-k) can expressed in term of (f/ij^., ■0k) from ^_k = 
r*;; and *L^ = *Jr with r = CTi. Then we get that 
t/(k) has an additional constraint C/(k) = — r[/-^(k)r. 
Under the constraint, we divide 4 Emit 2x2 matrices 
Ma = (Ijf) into two classes : in one class, the matrix 
M3 = (T3 commutes with F, we call it " even matrix^^; the 
others Mq = 1, Mi = ci and A/2 = 172 anti-commute 
with F, we call it "odd matrices". Then one can ex- 
pand C/(k) as C/(k) = 2u„(k)Afc<, a = 0,1,2,3. We 

a 

can prove that u^{\i) ~ U3(— k) and Ua{\i) — — Uq(— k), 
a = (0, 1, 2). Thus Ua(k) ( a = (0, 1, 2) ) are odd func- 
tions of kx, ky and are fixed to be zero at momenta (0, 0), 
(0,7r), (7r,0), (tt.tt) : Ua(k = 0) = 0. Here k = de- 
note four special points in momentum space (0, 0) , (0, tt), 
(7r,0), (7r,7r) . 



In general, ?7(k) is written into 

TTC\^\ _ Z' £k Ai,k + iA2.k 

^W-1 Ai^k-*A2,k -Ek 



where e^ is the even function of k^ and ky and Aik, A2,k 
are the odd functions of k^ 
one has 



and ky. At the points k > 0, 



H{k > 0)r 



E 

k>0 



[alpl 



e{k) 




-e{k) 



"k 



where ak>o and /?k>o == Q^lk are (diagonahzed) quasi- 
particle operators. Here k > means that fc^ > or 
ky = 0, k^ > 0. At the four points (0,0), (0,7r), (tt, 0), 
(tt, tt), the "mean- field" Hamiltonian is diagonalizcd into 



H{k = 0)r 



k=0 



e{k)ipliPi 



kVk- 



The ground state is denoted by a projected p-wave BCS 
type wave-function as 

1.9) = \^) = i^ni"kr/'exp(i^.gk^tV^k)|0). (2) 



where gk = v^/uk and 



2 \^-^ _Bk 



1 



Bk 



Here ^k = V4 + (^i.k)' + (A2,k)2 is the 
energy dispersion. P is the projected operator which 
guarantees a constraint to the Hilbert space with single 
fcrmion on each site. For the fermionic vacuum we have 
V'k|0) = and all the negative energy levels are filled by 
the fermions. Under these conditions, the ground state 
of Eq.([2]) is not a p-wave superconducting state, instead 
it is an insulator of topological order! 

In the following part we will use the projected p-wave 
SC wave functions to learn the nature of the topological 
orders, including the classification of topological orders 
for spin models with translation invariance, the degener- 
acy of the ground states and the prediction of new types 
of topological orders. 

Classification of topological orders - Firstly, we clas- 
sify the topological orders on the lattice by the projected 
p-wave SC wave-functions. Because there is no local or- 
der parameter to distinguish different topological phases, 
we classify the topological phases by the quantum phase 
transitions between them. Since the quantum states with 
topological order are protected by the finite energy gap 
for excitations and stable against any local perturbations, 
to break down a topological phase, one needs to close 
the energy gap for the excitations [9|, llO|. That is the 
unmovable gapless excitations indicate a quantum phase 
transition between different topological orders. 

For the topological orders described by the projected 
p-wave SC wave-functions, we find that e(k = 0) = 
at one or more points of k = define the quantum 




FIG. 1; Classification of topological orders by a p-wave SC 
wave function for Jz > 0- 




FIG. 2: Classification of topological orders by a p-wave SC 
wave function for J, < 0. 



topological transitions between different topological or- 
ders. For example, the topological order describe by the 
projected p-wave SC wave-function with e(k = 0) > 
and that by the projected p-wave SC wave-function with 
e(k = 0) < cannot be the same one. The en- 
ergy gap must vanish at k = 0, corresponding to the 
gapless critical states separating two topological phases. 
Then the sighs of e{k = 0) become an additional struc- 
ture beyond projected spin group to distinguish topolog- 
ical phases. There are totally 16 different cases for the 
sighs of £(k ~ 0) which characterize 16 different universal 
classes of topological orders. We use the symbols "/ii^A^" 
to denote them. /i(or i/. A, ^ ) = 1 denotes the case for 
£(k = (0, 0)) < (or e(k = (tt, 0)) < 0, e(k = (0, tt)) < 0, 
e(k = (tt, tt)) < 0) and /i(or i/. A, ^ ) = denotes the case 
e(k = (0, 0)) > (or e(k == (tt, 0)) < 0, e(k = (0, tt)) < 0, 
e(k==(^,7r))<0). 

Fermionic wave-functions for the topological orders - 
let's give one example for each class of topological order 
by the p-wave SC wave- functions. A general mean field 



p-wavc SC wave-function is written as 



f/(k) 



Ai k - iA2,k -£k 



S\^ — U 2; fJ X COS Kj; tjy COS f^y ^ 



Al,9 



J sin kx + J sin fe„, A2 k = 



(3) 



with J ^ 0. The classification of the topological orders for 
above wave-function is shown in the fig.l and fig. 2. From 
the two figures, one can sec that the topological orders 
of class 1001 and class 0110 cannot be denoted by above 
wave functions. For 1001, instead, a given ansatz is Sk = 
—J cos{kx + ky), Ai^k = J'sinfca;, A2,k = —J sin{kx + ky). 
For 0110, a given ansatz is Ck = J cos{kx + ky), Ai.k = 
J'sinfca;, A2,k = Jsin(fca; -|- ky). 

Topological degeneracies - We calculate the topolog- 
ical degeneracies for different classes of topological or- 
ders by using the p-wave SC wave-functions. Now 
we assume (u|7'"\ r;^™'"^) = ((_)™«xfe)(_)n.„fe)^^^.^ 
(^^yns^i'ij) f^_yisy{ij}^_-^ ^Q note four degenerate ground 
states for the topological orders on a torus. Take m,n — 
0,1 as an example. Sj.^y{ij) have values or 1, with 
Sx.y{ij) = 1 if the link ij crosses the x oi y line and 
Sx,y{ij) = otherwise. However, some ansatzs may be 
not permitted by the projections P^ (or P°). The ground 

state |\l/mean ' '^ } with evcn (or odd) number parti- 
cles is only allowed under projection P^ (or P°), 



For the mean-field ansatz defined above, there always 
exists a negative energy level for fermions at the point 
k > 0. For the energy of /3 band — £(k) is negative, 
£(k) > 0, /3 band is occupied and a band is empty, we 
have A'k>o = 0. On the other hand, for the energy of a 
band £(k) is negative, e(k) < 0, a band is occupied and 
f3 band is empty, we have A^k>o = 2. As a result, the 
total number of ip fermions on every point of k > 0, 
Y^ Nit., is even. Then (—1)^ is reduced into (—1)^" 

k>0 

with N = J2 Nk " 



Nf). Here N, 







J2 V'k'^k is the total 

k>0 k=0 

fermion number on the four special points k = 0. 



So to determine if the mean-field ground state con- 
tain even or odd number of "0 fermions, we only need 
to examine the occupation on the four special points: 
(0, 0), (0, tt), (tt, 0), (tt, tt) and calculate the total fermion 
number on them. The projection operators now are re- 



duced into PS 



i+(--i)^ 



and P? 



i-(-i)"» 



For a 



2 " r 2 

given point k = 0, the energy for 'ip is £(k = 0). Then if 

e(k = 0) < 0, one ^p particle occupies energy level at the 
point (0,0) (or (0,7r), (7r,0), (7r,7r)); if £r(k = 0) > 0, zero 
ip particle occupies the energy level at (0,0) (or (0, tt), 
(tt, 0), (tt, 7r)). Finally we obtain the projection opera- 
tors as 



1±(-1)^« _ 1±(-1)'^ 



E C(T(k)) 



Vp . \ — P'°>o|vl/ 



(«<r 



^.r 



j — I ^mcan 



Y (4) 



To obtain the topological degeneracy, one needs to deal 
with the projection operators -P^ and P° carefully. Let's 
simplify the projection operators firstly. We note that, 
for k > 0, 



^k>0 — Q^k>0'^k>0 - /3k>o/5k>0 + 1- 



(5) 



where the function C,{x) denotes C,(x) = 1, for x < 
and C {x) = 0, for x > 0. One can check the occupation 
number at k = easily to know if a given ansatz is valid. 
Using this method, we obtain the topological degeneracy 
on different lattices for the 16 topological orders. The 
results are given in the following tables in detail : 



nil 1110 1101 1011 0111 1100 1010 1001 OllO OOll OlOl lOOO OlOO OOIO 0001 0000 
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From the results, we can see that for the 16 topological 
orders, there are 8 Z2 topological orders ( 1111, 1100, 

1010, 1001, 0101, 0011, 0110, 0000) and 8 non-Abelian 
topological orders (1000, 0100, 0010, 0001, 1110, 1101, 

1011, 0111). 



The Z2 topological orders are described by Z2 projec- 
tive symmetry groups. Such Z2 topological orders always 
have 4 degenerate ground states on an even-by-even lat- 
tice with periodic boundary condition. And there always 
exists three types of quasiparticles: Z2 charge, Z2 vortex. 



and fcrmions, respectively [l]|. And the fcrmions can be 
regarded as bound states of a Z2 charge and a Z2 vortex. 

On the other hand, we have two classes of 8 types 
the non-Abelian topological orders [13 : on an even- by- 
even lattice the topological degeneracy is always 3. How- 
ever, on an odd-by-odd lattice, the degeneracy is 3 for 
first class 4 non-Abelian topological orders and 1 for the 
other class 4 topological orders. All the 8 non-Abelian 
topological orders have the two types of quasiparticles: 
non-Abelian anyons and fermions. 

Our results are consistent to those by Read and Green 
[l3|. They have pointed out there are two kinds of 
Px + ipy-yfsve fermion paired state : one is in the weak 
pairing phase as an example of topological order with 
non-Abelian anyon [13|, the other in the strong pairing 
phase as a topological order with Abelian anyon. How- 
ever, our results show more rich physics properties for 
the Px + ipy-w&Ye fermion paired states on lattice beyond 
the picture of weak and strong pairing phases. 

Unification of the topological orders for exact solved 
spin models - Finally we can identity the classes of topo- 
logical orders for the exact solved spin models by our 
classification and show their topological properties. 

For the Kitaev model on a hexagonal lattice [J], the 
original spin model is mapped to a model of p -wave BCS 
model on a square lattice by a Jordan- Wigner transfor- 
mation with [5| 

£k = |Jz| - |J2:|c0sfc:i, - |Jy|cOSfcy, 

Ai,5 = \Jx\sinkx + \Jy\sinky, A2^k^0. (6) 

From the relationship between | J^,] , | Jy|, | J^j , we found 
that the phase \Jz\ > \Jx\ + \Jy\ as 0000 type Z2 topo- 
logical order; \Jy\ > \Jx\ + \Jz\ as 1010 type Z2 topo- 
logical order; \Jx\ > \Jy\ + \Jz\ as 1100 type Z2 topo- 
logical order. The ground state degenerate is always 4 
on an even- by-even lattice (on a torus). However, on 
other lattices (even- by-odd, odd-by-even, odd-by-odd), 
the ground state degenerate is different : for 0000 type, 
it is 4; for 1010; it is 4 on an odd-by-even lattice, but 2 
on an even-by-odd and odd-by-odd lattice; for 1100; it 
is 4 on an even-by-odd lattice, but 2 on an odd-by-cvcn 
and odd-by-odd lattice. 

For the Kitaev model on a honeycomb lattice with min- 
imal three- and four-spin terms and a T-symmetry break- 
ing external magnetic field, the corresponding p- wave SC 
state is given as |6[ [Zl |8| 



3, 3, 1 for even-by-even, even-by-odd, odd-by-cvcn, odd- 
by-odd lattices, respectively. 

For the Wen's plaquette model, a "mean-field" ansatz 
p-wave BCS states bccomes[2,l3J 



H„ 



E (^KV'j., + ^Ivi/^h + ''■^) (8) 



{ij) 



where I,J = 1,2. The ground state for g < is a 



topological state with —rji^i^x 



^i,i-\-x 



1 4- cr and 



Vi,i+y 



^i,i~^y 



1 — c . The quantum wave- function 



can be regarded as two-flavor spinless projected p-wave 
SC. Then the fermion number on each lattice and the 
total fermion number must be even. Using the same 
method, one can also calculate the ground state degen- 
eracy which is 4, 2, 2, 2 for even-by-even, even-by-odd, 
odd-by-even, odd-by-odd lattices, respectively. It is the 
0110 type of the topological order. 

In addition, it is known that the topological orders of 
the toric-code model and the Wen-plaquette model are 
the the same ones[l|. So one can use the same p-wave 
SC wave-function in Eq.(8) to describe the toric code 
model. 

In conclusion, we develop a sysmatical theory for 
topological orders by the (projected) p-wave SC wave- 
functions. Based on the theory we unify different topo- 
logical ordered states for the exact solved spin models 
and obtain their topological properties. From our classi- 
fication in this paper, the topological orders in spin mod- 
els are the same one and at least five different classes of 
topological orders have been found. That is one can- 
not change the topological state for the Wen-plaquette 
model into that of the Kitaev model without a quantum 
phase transition. And we can predict that at the quan- 
tum phase transition between them the massless excita- 
tions should have zero energy at k = 0. In the end, we 
can also give a prediction that there are totally 11 un- 
known classes topological orders, including four classes 
of Z2 topological order ( 1111, 1001, 0011, 0101 ) and 
7 classes of non-Abelian topological orders. All the new 
classes of topological orders need to be explored in the 
future. 
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Ek 



\Jz 



Jx 



COSfcr 



J^ 



COSfcy, 



Ai.k = Aix sin kx + Aiy sin ky , 
A2,k = A2X sin kx + A2y sin ky 



with \JJ < \Jx\ 



^V 1 ^ z ^ ^x 



^V 1 ^ z 



(7) 



> 



— I >/x I + 1 >/j/ 1 ■ It is described by the 1000 type Non-Abelian 
topological order[4|,|l4|. The topological degeneracy is 3, 
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